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Abstract 

A  thermo-electro-hydro-dynamic  model  is  developed  to  examine  the  effects  of  fluid  slip  on  electrokinetic  energy  conversion  in  nanofluidic 
channels.  The  optimum  electrokinetic  devices  performance  is  dependent  on  a  figure  of  merit  defined  in  terms  of  hydrodynamic  conductance, 
streaming  effects  and  electrical  conductance  of  slip  nanochannels.  These  phenomenological  coefficients  are  all  enhanced  by  the  fluid  slip  as 
compared  to  those  without  slip  effects.  The  net  result  is  an  increased  figure  of  merit  and  thus  an  enhanced  electrokinetic  devices  performance, 
particularly  in  nanochannels  with  a  high  ratio  of  slip  length  to  channel  height. 

©  2007  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Recent  research  in  the  field  of  electrokinetic  energy  conver¬ 
sion  has  focused  on  increasing  the  efficiency  of  electrokinetic 
generation  [1-3]  or  pumping  [4-6]  in  nanofluidic  channels. 
Though  a  maximum  efficiency  of  15%  has  been  theoretically 
predicted  [4,7],  the  experimentally  measured  efficiencies  of 
electrokinetic  energy  conversion  have  typically  been  less  than 
1%  [8-11],  leaving  a  large  room  for  improvement  of  elec¬ 
trokinetic  devices  performance.  More  recently,  Pennathur  et 
al.  [12]  proposed  the  use  of  slip  nanochannels  [13]  for  elec¬ 
trokinetic  energy  conversion.  They  claimed  that  a  conversion 
efficiency  of  up  to  35%  is  achievable  in  lOOnm-diameter  cylin¬ 
drical  tubes  with  a  slip  length  of  6.5  nm  and  a  very  high  surface 
potential.  However,  the  real  effects  of  fluid  slip  on  the  flow  of 
liquid  and  ions  (i.e.,  electric  current)  and  thus  the  electroki¬ 
netic  energy  conversion  in  nanofluidic  channels  still  remain 
unclear.  This  communication  aims  to  address  this  issue  using 
a  thermo-electro-hydro-dynamic  model  developed  by  Xuan  and 
Li  [7]. 
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2.  Theoretical  formulation 

2.7.  Nomenclature 


b  slip  length 

B  non-dimensional  slip  length 

Cb  ionic  concentration 

e  charge  of  proton 

gi  defined  non-dimensional  functions  (/  =  1,  2,  3) 
G  hydrodynamic  conductance 

h  half  channel  height 

i  electric  current  density 

I  electric  current 

Boltzmann’s  constant 
K  non-dimensional  channel  height 

M  phenomenological  coefficient 

Na  Avogadro’s  number 

p  hydrodynamic  pressure 

Q  fluid  volume  flow  rate 

R  universal  gas  constant 

S  electrical  conductance 

T  liquid  temperature 

u  fluid  velocity 

y  transverse  coordinate 

zv  valence  of  electrolyte  ions 

Z  figure  of  merit 
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Greek  letters 

P  non-dimensional  property  of  fluid  (Levine  number) 

s  fluid  permittivity 

f  zeta  potential 

normalized  zeta  potential 
rj  energy  conversion  efficiency 

k  inverse  of  Debye  screening  length 

kb  molar  conductivity  of  fluid 

ti  fluid  viscosity 

0  electric  potential 

0  electric  double-layer  potential 

E  normalized  double-layer  potential 

double-layer  potential  at  channel  center 

2.2.  Thermodynamic  analysis 

Electrokinetic  energy  conversion  devices  are  thermodynam¬ 
ically  characterized  by  the  Onsager  reciprocal  relations  for 
volume  flow  rate  Q  and  electric  current  I  [14] 

Q  =  G(—Vp)  +  M(— V0)  (1) 

/  =  M(-Vp)  +  S(-V0)  (2) 


profile  for  an  axial  slip  flow 


h 2  (  y2  2 b\ 

u  =  rAl  ^  + 


+ 


l1 


d/  b  d  0 

1 - 1 - - — 


(-V0) 


_y=/z  y 


(5) 


where  /r  is  the  fluid  viscosity,  £  is  the  fluid  permittivity,  £  is  the 
zeta  potential  on  the  channel  wall,  0  is  the  electric  double-layer 
potential.  It  is  apparent  from  Eq.  (5)  that  the  fluid  slip  b  serves  to 
enhance  both  the  pressure  driven  (1st  term)  and  electro-osmotic 
(2nd  term)  flows. 

The  current  density,  i,  in  a  slit-type  nanochannel  flow  is  given 
by  [17,18] 

d2\Lr  ( Zve\k\ 

i  =  ~£  Jy2u  +  CbA<b  cosh  ( ~j—f  J  (_V<W  (6) 

where  Cb  is  the  ionic  concentration  of  the  bulk  electrolyte,  kb 
is  the  molar  conductivity  of  the  electrolyte,  zv  is  the  valence 
of  electrolyte  ions,  e  is  the  charge  of  proton,  k b  is  the  Boltz¬ 
mann’s  constant,  and  Lis  the  fluid  temperature.  The  double-layer 
potential  0  is  solved  from  the  Poisson-Boltzmann  equation  [16] 


where  Vp  and  V0  are  the  pressure  gradient  and  electric  potential 
gradient  across  the  channel,  G  is  the  hydrodynamic  conductance, 
M  characterizes  the  electro-osmotic  flow/streaming  current, 
and  S  is  the  electrical  conductance.  Using  these  two  phe¬ 
nomenological  equations,  Xuan  and  Li  [7]  demonstrated  that  the 
optimal  performance  of  electrokinetic  devices  is  solely  depen¬ 
dent  on  Z,  a  non-dimensional  “figure  of  merit”  [15]  defined 
as 


Z  = 


M2 

~GS 


(3) 


Moreover,  the  maximum  efficiencies  of  electrokinetic  gener¬ 
ators  and  pumps,  r\ max,  are  identical  and  given  by  [7] 


Z 

^max  -  2(2  _  z) 


(4) 


As  the  figure  of  merit  Z  is  less  than  1  due  to  the  irreversibility, 
OmsL x  would  not  be  able  to  reach  50%.  The  phenomenological 
coefficients  G,  M  and  S  are  functions  of  liquid  and/or  channel 
properties,  and  can  be  specified  using  the  electrokinetic  flow 
theory. 


2.3.  Electro-hydro-dynamic  analysis 


dV 

Ay2 


=  k2  sin/z(d/) 


(7) 


where  -Zye^lk^T  is  the  normalized  double-layer  potential, 
k  =  yj 2Zye2  N\Cb /sk#  T  is  the  inverse  of  the  Debye  screening 
length  with  Na  the  Avogadro’s  number.  The  analytical  solution 
of  E  is  expressed  in  terms  of  a  Jacobian  elliptical  function  [19] 


ip  =  ip0  +  2  In 


JacCD  (  — e‘^l/2 

e2*°)} 

V  2 

J1 

(8) 


where  is  the  double-layer  potential  at  the  channel  center  and 
can  be  specified  using  the  known  zeta  potential  on  the  channel 
wall. 

Integrating  Eqs.  (5)  and  (6)  over  the  channel  cross-section 
and  then  comparing  the  resultant  equations  with  Eqs.  (1)  and  (2) 
yield  the  expressions  of  the  phenomenological  coefficients 


2hs£ 

M  = - - 

d 


(>-«+£«) 


4  s2CK2 
fxht;*2 


g3  -  cosh(tf'o)  + 


b_g\ 

h  2 


(9) 

(10) 

(ID 


Consider  a  slit  nanochannel  of  height  2 h,  width  w  and  length  / 
where  w^>2h  and  l  2 h,  so  that  the  end  effects  can  be  neglected 
and  the  channel  can  be  treated  as  infinite  parallel  plates.  To 
account  for  the  fluid  slip,  we  utilize  Navier’s  slip  boundary  con¬ 
dition,  u(h)  =  —b  du/dy\y=h,  where  u(h)  is  the  slip  velocity  at  the 
wall,  b  is  the  slip  length,  and  y  is  the  transverse  coordinate  orig¬ 
inating  from  the  channel  center.  Applying  this  slip  condition  to 
a  general  electrokinetic  flow  [16,17],  we  arrive  at  the  velocity 


where  8i=JoWOd(y/h),  g2  = 

V2[cosh(f*)  —  coshfd'o)]  and  g2  =  f  cosh(d/)d(y//!) 
are  all  positive  functions  [20],  -Zyet.lk^T  is  the  normalized 
zeta  potential,  K-ich  is  the  non-dimensional  channel  height, 
/3  =  kbp,/sRT  (previously  termed  Levine  number  [5])  is  the 
non-dimensional  property  of  the  bulk  electrolyte  with  R  the 
universal  gas  constant.  One  can  see  that  fluid  slip  enhances  all 
three  phenomenological  coefficients. 
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Substituting  Eqs.  (9)— (1 1)  into  Eq.  (3)  leads  to  the  figure  of 
merit  Z  in  slip  nanochannels 

z  =  3[(i-giK7^  +  gg2]2 

2(1  +  3B)[(1  +  0/4)g3  -  cosh(«B0)  +  Bgj/2] 

where  B  =  b/h  is  the  non-dimensional  slip  length.  Therefore, 
Z  is  a  function  of  four  non-dimensional  parameters,  viz.  /3, 
B ,  K  and  f*.  Apparently,  a  smaller  results  in  a  higher  Z. 
The  effect  of  B  on  Z,  however,  cannot  be  readily  inferred 
from  Eq.  (12)  because  it  causes  an  increase  in  each  of  the 
terms  in  both  the  numerator  and  the  denominator.  In  addition, 
Z  seems  to  favor  small  K  while  large  f*.  All  these  aspects 
will  be  analyzed  and  discussed  below,  where  1  mM  KC1  aque¬ 
ous  solution  (i.e.,  Cb  =  1  mM)  is  used  as  the  electrolyte.  The 
fluid  properties  include  viscosity  /x  =  0.9  x  10-3  kgm-1  s_1, 
permittivity  s  =  79  x  8.854  x  10-12  CV-1  m-1,  and  molar  con¬ 
ductivity  A,b  =  0.0144  S  m2  mol-1  at  temperature  T  =  298  K  [1 1]. 
The  Levine  number  /3  is  thus  7.47.  A  MATLAB  program  was 
written  to  implement  the  calculations.  An  iterative  method  was 
first  employed  to  determine  the  double-layer  potential  from 
Eq.  (8).  Then,  a  direct  numerical  integration  approach  was 
applied  to  compute  the  phenomenological  coefficients  and  figure 
of  merit  in  slip  nanochannels. 

3.  Results  and  discussion 

Fig.  1  displays  the  ratios  of  the  phenomenological  coeffi¬ 
cients  G,  M  and  S  in  nanochannels  with  a  slip  length  of  b  =  5  nm 
to  those  without  slip  (i.e.,  b  =  0).  As  expected,  each  coefficient  is 
enhanced  by  the  effects  of  fluid  slip,  where  the  largest  enhance¬ 
ments  are  obtained  by  G  and  M  in  the  region  of  small  K  values 
(or  equivalently,  small  nanochannels  with  strong  double-layer 
overlapping)  while  S  remains  mostly  constant.  The  net  effects 
of  fluid  slip  are  still  to  increase  the  figure  of  merit  Z  as  illustrated 
in  Fig.  2.  For  example,  the  maximum  Z  value  in  slip  nanochan¬ 
nels  (top)  is  0.83  as  compared  to  0.24  in  no-slip  nanochannels 
(bottom).  Moreover,  fluid  slip  squeezes  the  region  within  which 
Z  approaches  its  extreme  to  the  one  with  a  higher  magnitude  of 
£*  while  a  lower  value  of  K. 


A 


Fig.  1 .  Ratios  of  the  phenomenological  coefficients  G,  M  and  S  in  nanochannels 
with  a  slip  length  b  =  5  nm  to  those  without  slip  (i.e.,  b  =  0).  The  normalized  zeta 
potential  is  £*  =  4. 


1  2  3  4  5  6 
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Fig.  2.  Contours  of  the  figure  of  merit  Z  as  a  function  of  the  normalized  zeta 
potential  £  *  and  the  non-dimensional  channel  height  K  when  fluid  slip  (b  =  5  nm) 
is  (top)  considered  and  (bottom)  ignored. 

Fig.  3  compares  the  maximum  efficiency,  ??max  in  Eq.  (4),  of 
electrokinetic  energy  conversion  with  (solid  lines)  and  without 
(dashed  lines)  consideration  of  fluid  slip  (b  =  5  nm)  at  differ¬ 
ent  £*  values.  While  ?7max  increases  for  all  values  of  £*  and 
K  due  to  the  effects  of  fluid  slip,  it  can  be  clearly  seen  that 
narrow  channels  with  high  zeta  potentials  benefit  the  most.  A 
maximum  efficiency  of  more  than  30%  is  obtained  in  the  slip 
nanochannel  with  AT~0.3  (corresponding  to  a  normalized  slip 
length,  #  =  0.26)  and  £*  =  6.  Such  high  conversion  efficiency 
will  make  the  electrokinetic  energy  conversion  devices  in  slip 
nanochannels  very  competitive  and  attractive.  As  a  comparison, 
the  maximum  efficiency  in  no-slip  nanochannels  is  only  7.2%  at 
K  ~  0.9,  which  is  less  than  a  quarter  of  that  in  slip  nanochannels. 


A 

Fig.  3.  Comparison  of  the  electrokinetic  conversion  efficiency  r]max  when  fluid 
slip  (b  =  5  nm)  is  considered  (solid  lines)  and  ignored  (dashed  lines). 


300 


C.  Davidson ,  X.  Xuan  /  Journal  of  Power  Sources  179  (2008)  297-300 


t  i  i  i  i  i 

0  2  4  6  8  10 

b  (nm) 


Fig.  4.  Electrokinetic  conversion  efficiency  ijmax  (normalized  by  that  with  zero 
slip  length)  as  a  function  of  the  non-dimensional  slip  length  B.  The  dimensional 
slip  length  b  is  also  shown  as  a  second  abscissa.  The  non-dimensional  channel 
height  is  fixed  to  K = 2. 

Moreover,  as  demonstrated  in  Fig.  2  (top),  an  even  higher  zeta 
potential  will  further  increase  the  Z  value  and  thus  ^max  in  Fig.  3. 

The  influence  of  slip  length  on  ^max  is  displayed  in  Fig.  4, 
where  r]mSLX  is  normalized  by  that  without  fluid  slip  and  K  is 
fixed  to  2.  The  efficiency  increases  with  the  non-dimensional  slip 
length  B  (and  thus  the  dimensional  slip  length  b  as  the  channel 
half  height  h  is  fixed;  b  is  also  shown  in  Fig.  4  as  a  second 
abscissa),  and  is  approximately  a  linear  function  of  B  especially 
when  £*  is  small.  This  phenomenon  may  be  explained  by  the 
fact  that  the  increase  in  M  due  to  the  effects  of  fluid  slip  is  almost 
cancelled  out  by  the  increases  in  G  and  S  (see  Fig.  1).  As  such, 
the  slip  effect  on  Z  is  roughly  equivalent  to  that  on  M  where  B 
appears  only  in  the  linear  term  [see  Eq.  (12)]. 

4.  Conclusions 

We  have  examined  the  effects  of  fluid  slip  on  electrokinetic 
energy  conversion  in  nanochannels  using  an  analytical  model 
based  on  Onsager’s  reciprocal  relations  and  electrokinetic  flow 
theory.  Fluid  slip  is  found  to  enhance  the  phenomenological 


coefficients  and  thus  the  electrokinetic  figure  of  merit.  The  result 
is  an  approximately  linear  increase  in  the  electrokinetic  devices 
performance  with  respect  to  the  slip  length.  Our  results  indicate 
that  nanochannels  with  a  larger  ratio  of  slip  length  to  channel 
height  will  benefit  more  from  the  fluid  slip  effects,  particularly 
when  the  wall  surface  potential  is  high. 
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